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Uncertainty  Analysis  of  the  NIST  Nitrogen  Flow  Facility 

Jennifer  L.  Scott  and  Michael  A.  Lewis 

National  Institute  of  Standards  and  Technology 

Chemical  Science  and  Technology  Laboratory 

Boulder,  Colorado  80303 

An  uncertainty  analysis  of  the  nitrogen  flow  facility  at  the  National 
Institute  of  Standards  and  Technology  was  performed.  This  facility  functions 
as  a  cryogenic  flow  calibration  laboratory  and  as  an  applied  research 
laboratory  for  high  pressure  nitrogen  gas  flow  measurement.  This  report 
includes  the  analysis  of  uncertainty  in  liquid  nitrogen  mass  flow,  gaseous 
nitrogen  mass  flow  determined  by  a  system  turbine  meter,  instrumentation, 
and  the  uncertainty  in  orifice  meter  discharge  coefficient  calculation  using  a 
defined  propagation  of  uncertainties  technique.  Uncertainties  determined  by 
statistical  means  are  distinguished  from  those  determined  by  other  means. 

Key  words:  buoyancy;  density;  discharge  coefficient;  mass;  propagation  of 
uncertainties;  orifice  meter;  sensitivity;  statistics;  turbine  meter;  Type  A 
uncertainty,  Type  B  uncertainty 

1.  Introduction 

This  report  provides  a  detailed  description  of  the  uncertainty  analysis  of  the  nitrogen 
flow  facility  at  the  National  Institute  of  Standards  and  Technology  (NIST)  in  Boulder, 
Colorado.  This  facility  is  used  in  two  distinct  ways.  We  calibrate  and/or  test  cryogenic 
flowmeters  using  a  mass  and  time  dynamic  weigh  system.  After  a  system  modification,  we 
can  conduct  applied  research  on  gas  flowmeters  using  the  same  weigh  system.  We  focus  the 
uncertainty  analysis  in  two  areas:  (1)  the  uncertainty  associated  with  the  determination  of 
mass  flow  rate,  and  (2)  the  uncertainty  in  measuring  orifice  meter  discharge  coefficients  in 
nitrogen  gas  flow. 

2.  Uncertainty  Evaluation 

We  use  the  guidelines  in  NIST  Technical  Note  1297  [1]  to  classify  the  types  of 
uncertainties  in  our  facility  and  to  determine  their  values.  We  distinguish  the  components 
of  uncertainty  determined  by  statistical  analysis  of  observed  data,  called  Type  A,  from  those 
determined  by  a  priori  information  (for  example,  manufacturers'  imcertainties  for  electronic 
equipment).  Type  B. 

Nearly  all  samples  used  to  estimate  Type  A  uncertainties  contain  more  than  30 
points.  Most  Type  B  uncertainties  are  those  associated  with  instrumentation  and 
thermophysical  properties,  of  the  working  fluids  (hquid  and  gaseous  nitrogen)  where  the 
number  of  degrees  of  freedom  is  considered  infinite.  Following  the  guideUne  in  Technical 


Note  1297,  we  assume  that  uncertainties  stated  by  manufacturers  are  based  on  a  rectangular 
distribution;  that  is,  the  performance  of  the  instrument  will  fall  within  the  manufacturers' 
uncertainty  100  percent  of  the  time.  To  estimate  a  standard  uncertainty  of  an  instrument, 
we  divide  the  manufacturer's  uncertainty  by  the  square  root  of  3. 

The  combination  of  all  uncertainties,  in  quadrature,  is  multiplied  by  a  coverage  factor 
of  2  in  order  to  estimate  an  approximate  95  percent  confidence  interval  for  total  system 
uncertainty.  All  sample  sizes  were  large  enough  that  a  coverage  factor  of  2  is  sufficient  to 
estimate  a  95  percent  confidence  interval. 

The  key  to  determining  the  uncertainty  in  our  system  centers  on  our  ability  to 
determine  the  uncertainty  in  crucial  measurements  of  temperature,  pressure,  frequency,  and 
mass.  Certain  electronic  instruments  are  common  to  many  of  these  measurements.  Our 
analyses  of  the  uncertainties  in  these  instruments,  as  well  as  the  uncertainties  in  density 
computation  and  dimensional  measurements,  are  included  in  the  Appendices.  The 
Appendices  also  contain  uncertainty  analyses  of  pressure  and  temperature  measurements. 
The  information  in  these  Appendices  lays  the  groundwork  for  determining  uncertainties  in 
measurements  involving  combinations  of  these  values. 

A  first-order  Taylor  series  expansion,  commonly  called  propagation  of  uncertainty, 
is  used  to  determine  the  uncertainties  in  a  quantity  (such  as,  discharge  coefficient)  that  is 
calculated  when  the  variables  in  the  calculation  have  associated  uncertainties.  Equations 
(1)  and  (2)  show  the  general  form  used  in  evaluating  propagation  of  uncertainties, 

y  =  /  (J^p  x^y-.jc^)  ,  (1) 


«'(y)  =  E 


i=l 


6x.. 


n-l 


"'i^,)  -  2E  E  ^^  u  (X,  X)  ,  (2) 


where  y  represents  a  quantity  that  is  calculated  as  a  function  of  n  other  quantities,  Xi,...Xn, 
and  Uc(y)  is  the  combined  standard  uncertainty  of  the  calculated  quantity.  The  third  term 
in  Equation  (2)  represents  the  covariance  between  variables. 

To  determine  uncertainties  in  many  components  of  the  flow  system  we  combine 
factors  as  though  their  relationship  was  multiplicative;  that  is,  we  sum  the  squares  of  the 
individual  uncertainties.  Whenever  the  functional  relationship  can  be  determined,  we  use 
the  propagation  of  uncertainty.  We  do  not  include  covariances  because  in  no  case  did  we 
determine  that  the  uncertainty  of  one  factor  was  a  function  of  the  uncertainty  in  another 
factor.  The  following  list  describes  notations  used  in  this  report. 


Notation: 

Type  A  Uncertainty  determined  by  statistical  means 

Type  B  Uncertainty  determined  by  other  means 

n  Sample  size 

k    A  coverage  factor  that  is  multiplied  by  the  combined  standard 
uncertainty  to  achieve  a  95  percent  confidence  interval  estimate 
for  total  uncertainty.  This  product  is  called  the  expanded 
uncertainty,  and  in  this  analysis,  k  equals  2. 

quadrature    (a^  +  b^  +  ....)^,  root  sum  of  the  squares 

u(Xj)    Estimate  of  the  standard  uncertainty  (la)  of  the  variable  Xj 

residual  standard    Standard  deviation  of  the  mean  of  the  difference  between  known 
deviation    values  and  those  predicted  by  a  model;  root  mean  squared  error 

standard  error    s//n,  standard  deviation  of  the  mean  divided  by  the  square  root  of 
of  the  mean    the  sample  size 

T    Temperature,  K 

V  Voltage,  V 
I    Current,  A 

p     Density,  kg/m^ 

CY  The  product  of  the  orifice  meter  discharge  coefficient  and  the 
expansion  factor.  This  product  is  referred  to  as  the  discharge 
coefficient  in  this  document. 

m  Mass,  kg 

t  Time,  s 

Ap  Differential  pressure;  pressure  drop  across  an  orifice  plate,  kPa 

static  pressure  Thermodynamic  or  line  pressure,  Pa 

D  Internal  pipe  diameter,  cm 

d  Orifice  bore  diameter,  cm 

Fa  Coefficient  of  thermal  expansion  for  stainless  steel 

/5  Beta  ratio,  ratio  of  the  orifice  bore  to  the  pipe  diameter  (d/D) 

test  point    This  phrase  describes  one  mass  measurement  by  either  the  weigh 
tank  or  a  turbine  meter.   It  may  include  any  or  all  of  the  following 
additional  measurements:   elapsed  time,  pulses,  temperatures,  and 
pressures. 

Y  Expansion  factor 


3.   Liquid  Nitrogen  Flow  Facility 
3.1  Facility  Description 

The  cryogenic  liquid  nitrogen  facility,  built  by  the  National  Bureau  of  Standards 
(NBS)  and  the  Compressed  Gas  Association,  was  commissioned  in  1968.  A  detailed 
description  of  the  early  capabilities  and  design  information  is  given  in  NBS  Report  9749  [2]. 
A  provisional  accuracy  statement,  NBS  Technical  Note  606  [3],  was  published  in  1971 
regarding  the  Hquid  nitrogen  system.  The  present  report  describes  slightly  different 
techniques  to  evaluate  uncertainty  but  relies  on  some  of  the  results  included  in  [3].  In 
addition,  changes  in  instrumentation  and  data  acquisition  since  the  pubUcation  of  Technical 
Note  606  require  a  new  uncertainty  evaluation.  This  section  describes  the  hardware  of  the 
hquid  facihty  and  the  measurements  required  for  flowmeter  caUbrations. 

Figure  1  is  a  schematic  of  the  liquid  flow  facility.  Liquid  nitrogen  is  the  process  fluid 
which  is  circulated  throughout  the  closed  loop  by  a  variable  speed  boost  pump.  The  liquid 
flows  into  a  subcooler  where  thermal  energy  due  to  pumping  and  ambient  heat  leak  is 
removed.  The  temperature  of  the  liquid  nitrogen  in  the  flow  loop  can  be  changed  or 
controlled  by  adjusting  the  liquid  level  and  the  vapor  pressure  in  the  subcooler  tank.  Some 
of  the  test  fluid  also  can  be  diverted  around  the  subcooler,  if  necessary. 

The  path  of  the  hquid  nitrogen  flow  is  shown  in  Figure  1.  After  leaving  the 
centrifugal  pump,  the  hquid  nitrogen  passes  through  the  subcooler  and/or  bypass,  through 
a  vacuum-jacketed  loop  to  the  test  section,  and  into  a  diffuser  which  removes  the  vertical 
component  of  velocity  at  the  bottom  of  a  0.378  m^  (100  gal)  aluminum  weigh  tank.  The 
hquid  flows  through  a  valve  at  the  bottom  of  the  weigh  tank  and  into  a  stainless  steel 
pressure  vessel  with  a  0.443  m^  (117  gal)  usable  capacity.  When  the  temperatures  and 
pressure  in  the  system  have  reached  a  steady  state  condition,  the  weigh  tank  valve  is  closed 
and  Hquid  nitrogen  accumulates  in  the  weigh  tank. 

Once  the  liquid  reaches  a  preset  level  in  the  weigh  tank,  a  test  point  begins.  A  timer 
is  initiated  and  a  computer  begins  storing  digitized  information  about  the  system  and  any 
flow  measurement  device  installed.  For  example,  if  a  flowmeter  with  a  frequency  output 
is  installed  in  the  test  section,  a  counter  begins  totalizing  the  meter  pulses  while  the 
computer  records  digitized  outputs  from  pressure  transducers  and  thermometers. 


A  second,  preset  liquid  nitrogen  level  marks  the  end  of  the  test  point.  Data  acquisition 
stops  and  information  from  the  timer  and  counters  is  sent  to  the  main  computer.  A  load 
cell  is  used  to  measure  the  mass  of  the  liquid  accumulated  in  the  tank.  The  load  cell  output 
is  recorded  at  the  beginning  and  end  of  the  test  point  so  the  difference  in  voltage  is  a 
measure  of  the  accumulated  mass.  The  load  cell  is  calibrated  by  suspending  NIST- 
calibrated  weights  on  the  load  cell  before  flowmeters  are  tested. 

Helium  is  used  as  an  over-pressurant  in  the  catch  and  weigh  tanks  to  control  the 
system  pressure.  The  system  is  always  over-pressurized  to  prevent  boiling  in  the  liquid 
nitrogen.  The  nitrogen  is  always  subcooled  by  10  to  15  K.  Helium  absorption  in  the  liquid 
nitrogen  is  expected  not  to  exceed  0.5  mole  percent  according  to  data  by  DeVaney,  Dalton, 
and  Meeks  [4].  No  test  points  are  taken  if  evidence  of  bubbles  is  detected  visually  through 
a  sapphire  window  located  near  the  meter  being  tested. 

The  continuous  flow  loop  allows  for  the  establishment  of  system  temperature, 
pressure,  and  flow  rate  with  the  capability  to  maintain  steady  conditions  for  substantial  time 
periods.   Typical  calibration  parameters  are: 

Flow  rate:  0.95  to  9.5  kg/s  (2  to  20  Ib/s) 
Pressure:   0.4  to  0.76  MPa  (60  to  110  psia) 
Temperature:    80  to  90  K  (144  to  162  °R) 

3.2   Uncertainty  in  Liquid  Nitrogen  Mass  Flow  Measurement 

The  heart  of  our  flow  measurement  facility  is  the  weigh  tank  system.  When  we 
calibrate  cryogenic  mass  flowmeters,  the  uncertainty  in  the  calibration  data  is  equal  to  the 
uncertainty  in  our  mass  flow  measurement.  The  uncertainty  in  a  calibration  of  a  cryogenic 
volume  flowmeter  depends  on  our  ability  to  measure  mass  flow  and  liquid  test-section 
density.  When  we  calculate  the  uncertainty  in  the  discharge  coefficient  of  an  orifice 
flowmeter,  the  uncertainty  in  mass  flow  is  an  integral  portion  of  the  propagation  of 
uncertainties  associated  with  this  calculated  quantity.  Consequently,  we  have  included  a 
detailed  description  of  our  calculation  of  uncertainty  in  mass  measurement.  The 
components  of  mass  uncertainty  include  those  from  the  load  cell,  the  buoyancy  corrections, 
time  measurement,  and  changing  mass  between  the  test  section  and  the  weigh  tank.  Percent 
uncertainties  are  based  on  a  181.4  kg  (400  lb)  test  draft. 

3.2.1   Load  Cell  Sensitivity 

The  mass  of  the  liquid  is  measured  by  a  load  cell.  The  load  cell  is  calibrated  by 
applying  a  series  of  certified  weights  on  the  load  cell  and  measuring  the  output.  At  every 
calibration  point  we  measure  both  the  excitation  voltage  (E^^^)  to,  and  the  output  voltage 
(Eq^)  from,  the  load  cell  before  and  after  the  weight  is  added.  Also,  we  measure  the 
pressure  of  the  ullage  surrounding  the  load  cell.  The  sensitivity  of  the  load  cell  at  that  point 
is  calculated  as  a  ratio  of  the  change  in  weight  to  the  difference  in  the  voltage  ratios  before 
and  after  the  weight  is  added.   Equation  (3)  illustrates  this  calculation: 


Sensitivity  =  <^'^"«"  '"  '^ 
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E  is  the  symbol  for  measured  voltage.  The  denominator  will  be  called  aE  in  the  remainder 
of  this  discussion  even  though  it  is  a  dimensionless  quantity. 

A  sample  of  144  of  these  sensitivities  was  fitted  by  regression  versus  aE  and  pressure. 
This  fit  creates  a  sensitivity  'surface'  which  is  called  the  sensitivity  equation  in  the  following 
discussion.  Even  though  it  is  nearly  flat  in  the  aE  plane,  it  is  important  that  this  variable 
be  included  in  the  regression  fit.  We  calibrate  the  load  cell  at  3  discrete  weights,  but  we 
must  interpolate  between  these  values  during  nitrogen  flow  measurement.  The  total 
uncertainty  in  the  calculation  of  the  sensitivity  of  the  load  cell,  which  is  direct  correlation 
to  uncertainty  in  mass  measurement,  has  four  components:  weights,  voltage,  sensitivity 
equation,  and  pressure. 

3.2.1.1  Weights 

The  load  cell  calibration  weights  are  four  nominal  113.4  kg  (250  lb)  brass  cylinders. 
The  calibration  documents  and  uncertainty  information  are  given  in  Appendix  A. 

3.2.1.2  Voltage 

Two  digital  voltmeters  are  used  to  measure  the  voltages  associated  with  the  load  cell: 
one  measures  the  nearly  constant  excitation  voltage  to  the  bridge  of  the  load  cell  and  the 
other  measures  the  voltage  from  the  leg  of  the  bridge.  All  calculations  using  these  voltages 
are  in  the  form  of  aE,  explained  in  Section  3.2.1.  This  quantity  is  a  difference  between 
voltage  ratios.  Because  we  are  using  the  voltages  in  this  manner  and  not  relying  on  the 
absolute  values  of  voltage,  we  do  not  include  the  uncertainties  associated  with  the  digital 
voltmeters.  Random  uncertainties  in  the  instruments  will  be  reflected  in  the  data  points 
used  in  the  calibration  fit  for  the  sensitivity  function. 

3.2.1.3  Sensitivity  Equation 

The  uncertainty  associated  with  the  equation  used  in  the  computation  of  sensitivity 
is  the  ratio  of  residual  standard  deviation  of  the  fit  to  a  nominal  value  of  sensitivity  and  is 
0.050  percent.  The  nominal  value  chosen  was  the  mean  of  the  144  data  points  used  in  the 
fit.  The  sensitivity,  a  measure  of  the  load  cell  output  as  a  function  of  the  change  in  added 
weight,  is  nearly  constant.  Areas  of  uncertainty  associated  with  the  load  cell  are  difficult  or 
impossible  to  evaluate  in  situ.  These  components  include  randomness  in  the  digital 
voltmeters  (see  Section  3.3.2),  the  repeatablilty  of  the  load  cell,  and  the  static  uncertainty 
in  the  calibration.  The  uncertainty  associated  with  the  equation  includes  many  of  these 
components. 


The  form  of  the  sensitivity  equation  also  accounts  for  dynamic  uncertainties  in  the 
load  cell  signal.  During  load  cell  calibration,  we  suspend  three  discrete  weights  from  the 
load  cell.  During  a  data  point,  the  weight  on  the  load  cell  is  varying  dynamically  as  the 
weigh  tank  is  being  filled  with  liquid  nitrogen.  We  do  not  evaluate  the  possible  uncertainty 
in  the  load  cell  due  to  these  dynamic  conditions.  The  only  way  to  quantify  these  dynamic 
uncertainties  would  be  through  a  reference  meter;  however,  the  magnitude  of  the  dynamic 
uncertainty  would  be  orders  of  magnitude  smaller  that  the  uncertainty  associated  with  any 
reference  meter. 

3^.1.4  Pressure 

The  uncertainty  associated  with  the  pressure  measurement  contributes  to  the 
uncertainty  in  the  sensitivity  because  it  is  an  independent  variable  in  the  sensitivity  equation. 
The  uncertainty  evaluation  for  the  ullage  pressure  transducer,  P9,  is  included  in  Appendix 
E.  Using  the  propagation  of  uncertainties  technique  shown  in  Section  2,  we  assessed  the 
additional  uncertainty  in  sensitivity  due  to  pressure  at  10'^  percent,  essentially  zero.  Table 
1  lists  the  uncertainties  associated  with  the  load  cell  sensitivity. 

Table  1.   Load  Cell  Sensitivity  Uncertainty  (la) 


Uncertainty  in  Sens. 

Type  A 

TypeB 

Equation 

0.050% 

Weights 

0.002% 

Total  (in  quadrature) 

0.050% 

0.002% 

3.2.2   Buoyancy  Correction 


The  buoyancy  of  the  mass  of  liquid  nitrogen  collected  in  the  weigh  tank  must  be 
evaluated  for  every  mass  measurement.  There  are  three  components  to  the  correction  in 
mass  measurement  due  to  buoyancy:  (1)  the  buoyancy  of  the  Hquid  nitrogen  accumulated 
during  the  data  point,  (2)  the  buoyancy  of  the  immersed  diffuser  and  pipe  in  the  liquid,  and 
(3)  the  change  in  buoyancy  of  the  liquid  nitrogen  accumulated  before  the  data  point  begins. 
The  first  two  are  reflected  in  the  equation 


B.C.  =  K,^  .  p 
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where  B.C.  is  the  correction  due  to  buoyancy,  V,iq  is  the  volume  of  the  liquid  nitrogen 
accumulated  during  a  data  point,  Punage  is  the  density  of  the  catch  tank  ullage,  Vpip^  is  the 
volume  of  the  diffuser  and  pipe  immersed  in  the  liquid  nitrogen,  and  p^^  is  the  density  of 
the  liquid  nitrogen.  The  propagation  of  uncertainty  for  the  buoyancy  correction  equation 
follows: 
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The  first  term  on  the  right  side  of  the  equation  represents  the  uncertainty  in  the  buoyancy 
due  to  hquid  nitrogen  density  as  it  relates  to  Hquid  nitrogen  volume;  the  second,  due  to  the 
catch  tank  ullage  density;  and  the  third  due  to  liquid  density  as  it  relates  to  the  immersed 
pipe  and  diffuser.  Each  of  these  components  will  be  addressed  in  the  following  sections,  as 
will  the  effect  of  the  change  in  buoyancy  of  the  liquid  nitrogen  accumulated  before  the  data 
point  begins.  The  uncertainty  calculations  that  follow  represent  the  effect  that  these 
uncertainties  have  on  a  nominal  181.4  kg  mass. 

322.1  Volume  of  Liquid  in  Weigh  Tank 

Our  ability  to  determine  the  volume  of  the  accumulated  liquid  nitrogen  is  a  direct 
result  of  our  ability  to  evaluate  the  pressure,  temperature,  and,  thereby,  density  of  liquid 
nitrogen.  The  uncertainty  in  the  pressure  (P9  and  barometer)  and  temperature  (T13) 
measurement  in  the  weigh  tank  are  found  in  Appendices  E  and  F.  Appendix  B  outlines  our 
uncertainty  evaluation  for  liquid  nitrogen  density.  It  includes  the  uncertainty  in  MIPROPS, 
the  property  software  we  use  [5,6],  as  well  as  the  uncertainty  in  density  due  to  uncertainties 
in  temperature  and  pressure  measurements.  Using  a  38-point  data  sample,  we  have 
evaluated  the  effect  of  these  density  uncertainties  on  a  nominal  181.4  kg  mass.  The  results 
are  shown  in  Table  2. 


Table  2.  Uncertainties  (la)  in  Mass  Measurement  Caused  by  Uncertainties  in  Liquid 

Density 


Source  of  Uncertainty  in 
Liquid  Density 

Percent  change  in  a  nominal  18L4  kg  mass 

Type  A 

TypeB 

MIPROPS 

0.004% 

6p/dT  p   (c7=0.5  K) 

0.006% 

6(>/dF  T  (P9,Baro) 

0% 

0% 

Total  (in  quad) 

0% 

0.007% 

3.2.2.2   Density  of  Nitrogen-Helium  Ullage 

The  two  most  important  factors  in  determining  the  density  of  the  He-N2  mixture  are 
the  temperature  of  the  ullage  in  the  volume  surrounding  weigh  tank  and  the  model  used  to 
estimate  the  density  of  this  mixture. 

We  use  the  Redlich-Kwong  model,  as  modified  by  Chueh  and  Prausnitz,  for  the 
binary  equation-of-state  calculation  [7].  This  model  uses  the  catch  tank  ullage  temperature 
and  pressure,  and  the  liquid  nitrogen  temperature  as  input.  When  this  model  was  adapted 
to  our  system,  empirical  tests  by  NIST  personnel  led  to  an  estimate  of  3  to  5  percent 
uncertainty  in  the  ullage  mixture  density  calculation  using  this  model.  Using  a  38-point 
sample  of  archival  data,  we  determined  the  effect  of  the  buoyancy  correction  on  a  181.4  kg 
mass  measurement  if  the  ullage  density  deviated  by  4  percent.  This  value  is  shown  in  Table 
3. 

We  also  determined  the  sensitivity  of  this  model  to  uncertainties  in  temperature  and 
pressure  measurement.  Three  temperature  measurements  are  involved  in  computing  the 
ullage  density:  the  two  ullage  thermometers  (TS71  and  TS72)  and  the  liquid  nitrogen 
thermometer  (T13).  The  liquid  nitrogen  temperature  is  required  to  determine  vapor 
pressure.  The  uncertainties  associated  with  these  thermometers  are  shown  in  Appendk  F, 
and  those  for  pressure  are  shown  in  Appendix  E.  These  uncertainties  in  temperature  and 
pressure  affect  the  ullage  density  calculation  and,  therefore,  the  mass  measurement.  The 
resultant  uncertainties  on  a  181.4  kg  mass  are  shown  in  Table  3. 
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Table  3.  Uncertainties  (la)  in  Mass  Measurement  Caused  by  Uncertainties  in  Ullage 

Density 


Source  of  Uncertainty  in 
Ullage  Density 

Percent  change  in  a  nominal  181.4  kg        mass 

Type  A 

Type  B 

Redlich-Kwong   (4%) 

0.063% 

dp/dT  p   (a=lK) 
TS71,  TS72 

0.018% 

(5p/(5T  p   (a=0.5  K) 
T13 

0.006% 

6p/dF  T  (P9,Baro) 

0.0001% 

0.0002% 

Total  (in  quad) 

0.0001% 

0.066% 

3.2.2.3  Volume  of  Diffuser  and  Pipe 

The  mass  must  also  be  corrected  for  the  buoyancy  of  the  diffuser  and  pipe  that  are 
immersed  in  the  collected  liquid.  The  uncertainty  in  the  correction  is  a  function  of  the 
volume  of  the  diffuser  and  pipe,  and  the  density  of  the  liquid  in  the  weigh  tank.  We  use 
the  maximum  diffuser  and  pipe  volume  that  can  be  immersed  in  the  liquid  to  provide  the 
most  conservative  estimate  of  uncertainty  due  to  the  buoyancy  correction  of  the  pipe  and 
diffuser.  For  that  reason,  we  do  not  evaluate  the  effect  on  the  buoyancy  correction  due  to 
our  uncertainty  in  measuring  the  volume  of  the  pipe  and  diffuser.  We  evaluated  the  effect 
of  a  0.5  percent  uncertainty  in  density  to  the  buoyancy  correction  for  the  pipe  and  diffuser 
and  assessed  its  contribution  to  the  uncertainty  in  a  181.4  kg  mass.  The  result  was  0.0035 
percent.  Because  we  did  not  determine  this  value  by  statistical  means,  it  is  a  Type  B 
uncertainty. 

3.2.2.4  Liquid  Volume  Accumulated  before  Measurement 

As  explained  in  Section  3.1,  a  test  point  begins  by  closing  the  valve  in  the  bottom  of 
the  weigh  tank.  When  the  accumulated  mass  reaches  a  preset  level,  'level  A,'  data 
acquisition  begin.  At  the  end  of  the  point,  'level  B,'  the  mass  is  determined  via  the 
difference  in  load  cell  signals  between  'level  A'  and  'level  B.'  Our  analysis  shows  that  the 
buoyancy  of  the  liquid  accumulated  before  the  test  point  begins  may  change  because  the 
catch  tank  ullage  temperatures  surrounding  this  volume  may  change.  We  correct  for  this 
change  in  buoyancy,  but  due  to  the  uncertainties  (Type  B)  associated  with  thermometers 
T13,  TS71,  and  TS72,  we  estimate  an  added  0.01  percent  uncertainty  in  measuring  a  181.4 
kg  mass  due  to  this  effect. 

Table  4  lists  all  the  uncertainties  in  mass  measurement  due  to  uncertainties  in 
buoyancy  corrections. 
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Table  4.  Mass  Uncertainties  (la)  Caused  by  Buoyancy  Corrections 


Mass  Uncertainty  due  to 
Buoyancy  Corrections 

Type  A 

TypeB 

Source:  Liquid  Density 

0.012% 

Source:  Ullage  Density 

0.0001% 

0.066% 

Source:  Diffuser/Pipe  Vol 

0.0035% 

Source:  Volume  Before  Point 

0.01% 

Total  (in  quadrature) 

0.0001% 

0.068% 

3.2.3  Time 

In  most  cases,  mass  flow  rate  is  used  for  meter  calibration.  A  multi-function 
datalogger  is  used  to  measure  elapsed  time.  The  uncertainty  in  time  measurement  is  equal 
to  the  resolution,  0.001  s.  The  contribution  of  time  to  the  uncertainty  in  mass  flow  rate  is 
evaluated  at  a  nominal  time  of  100  s.  This  is  a  Type  B  uncertainty. 

Table  5.  Uncertainty  (lor)  in  Mass  Flow  Rate  Caused  by  Time 


Time  Uncertainty 

TypeB 

0.001  s  @  100  s 

0.001% 

32A  Volume  between  Test  Section  and  Weigh  Tank 

There  is  a  finite  volume  of  liquid  nitrogen  in  the  piping  between  the  meter  being 
cahbrated  and  the  weigh  tank.  A  measure  of  the  density  gradient  in  the  pipe  is  the 
difference  between  fluid  temperature  at  the  meter  end,  Point  A,  and  the  weigh  tank  end  of 
the  volume,  Point  B.  We  did  not  evaluate  changes  in  density  due  to  pressure  drop  because 
there  is  no  reason  to  believe  that  the  pressure  drop  changes  during  a  data  point.  If  the 
density  gradient  in  this  pipe  remains  constant  during  the  data  point,  we  assume  that  the 
mass  collected  in  the  weigh  tank  is  equivalent  to  the  mass  passing  through  the  meter. 

We  evaluated  changes  in  the  density  gradient  in  this  volume  for  eight  data  points. 
The  average  number  of  thermometer  readings  per  data  point  was  165.  We  plotted  the 
difference  between  the  temperature  at  the  beginning  of  the  pipe  and  the  end  of  the  pipe. 
For  thermally  stable  flows,  this  difference  is  constant  and  the  slope  through  a  regression  fit 
of  these  data  is  zero.  A  nonzero  slope  indicates  that  the  density  gradient  and  the  mass  in 
the  pipe  are  changing. 

We  found  that  the  average  change  in  the  temperature  difference  between  Points  A 
and  B  for  the  eight  points  evaluated  was  0.067  K  which  translates  to  a  0.05  percent  change 
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in  liquid  nitrogen  density  (Appendix  B).  The  density  of  liquid  nitrogen  at  nominal  values 
of  85  K  and  586  kPa  is  773.02  kg/m^.  A  0.05  percent  change  in  this  density  times  a  nominal 
volume  of  the  piping  (0.05  m  )  results  in  a  0.019  kg  mass  change  in  the  piping.  The 
resultant  uncertainty  in  a  181.4  kg  mass  measurement  is  0.011  percent. 

We  did  not  include  the  uncertainties  associated  with  the  thermometers  used  in  the 
evaluation.  They  are  included  in  other  areas  where  the  absolute  values  of  the  temperature 
are  needed.  In  this  case,  only  the  change  in  the  difference  between  temperature 
measurements  was  evaluated  and  not  the  differences  themselves.  The  random  variation  of 
the  thermometers  were  accounted  for  in  the  regression  fit  of  the  temperature  differences. 

The  following  table  exhibits  all  uncertainties  associated  with  our  measurement  of 
mass  and  mass-flow  rate. 

Table  6.   Total  Uncertainties  (la)  for  Mass  Measurement 


Source  of  Uncertainty: 

Type  A 

Type  B 

Combined 

Load  Cell  Sensitivity 

0.050% 

0.002% 

0.050% 

Buoyancy  Correction 

0.0001% 

0.068% 

0.068% 

Mass  between  Test  Section  and 
Weigh  Tank  (liquid  Nj) 

0.011% 

Total  for  Mass  Measurement 

0.051% 

0.068% 

0.085% 

Expanded  Uncertainty,  k=2 

0.170% 

Time 

0.001% 

0.001% 

Total  for  Mass  Flow  Rate 

0.051% 

0.068% 

0.085% 

Expanded  Uncertainty,  k=2 

0.170% 

3.3   Liquid  Volume  Flow  Rate  Measurement 

Our  abiHty  to  accurately  measure  liquid  nitrogen  volume  flow  rate  through  the  test 
section  depends  not  only  on  our  ability  to  measure  mass  flow  rate,  but  also  on  our 
determination  of  liquid  density.  The  transducers  used  to  evaluate  the  pressure  and 
temperature  of  the  Hquid  in  the  test  section  are  P7  and  TIO.  Their  uncertainties  are  listed 
in  Appendices  E  and  F.  The  uncertainties  associated  with  density  calculations  are  shown 
in  Appendix  B.  Table  7  indicates  the  total  uncertainty  in  determining  liquid  volume  flow 
rate.   Density  uncertainties  were  evaluated  for  liquid  nitrogen  at  85  K  and  620  kPa. 


13 


Table  7.  Uncertainty  (1<7)  in  Volumetric  Flow  Rate  Measurement 


Source  of  Uncertainty  in 
Liquid  Volume  Flow 

Nominal  Values:  T  =  85  K,  P  =  620  kPa, 
Mass  =  181.4  kg,  Time  =  100  seconds 

Type  A 

Type  B 

MIPROPS 

0.25% 

dp/dT\^ 

0.025% 

0.013% 

6p/d?  T  (P7,Baro) 

0.0002% 

0.0002% 

Uncert.  in  mass  flow  rate 

0.051% 

0.068% 

Total  (in  quad) 

0.057% 

0.259% 

Total  (Type  A+B,  in  quad),  k=2 

0.531% 

Without  MIPROPS  uncertainty,  k  =  2 

0.179% 

The  major  contributor  to  the  uncertainty  of  hquid  nitrogen  volumetric  flow 
measurement  is  the  uncertainty  in  the  liquid  density.  The  greatest  uncertainty  in  density 
comes  from  the  stated  uncertainty  in  the  equation  of  state  (MIPROPS).  When  the  objective 
of  using  the  uncertainty  statement  is  in  the  mediation  between  seller  and  buyer,  where  both 
have  accepted  the  same  values  for  density  as  determined  by  the  NIST  standard,  the 
uncertainty  in  density  due  to  the  equation  of  state  need  not  be  considered.  In  that  case,  the 
expanded  uncertainty  in  volumetric  flow  measurement  for  this  laboratory  (with  k=2)  is  0.18 
percent. 

4.   Gaseous  Nitrogen  Flow  Facility 

The  gaseous  nitrogen  flow  loop  has  been  used  almost  exclusively  for  applied  orifice 
meter  research.  Much  of  the  following  discussion  centers  around  our  determination  of 
uncertainty  in  orifice  meter  discharge  coefficient. 

4.1   Facility  Description 

A  schematic  of  the  gaseous  nitrogen  flow  facility  is  shown  in  Figure  2.  All  of  the 
components  of  the  liquid  system  are  used  but,  the  flow  is  diverted  as  it  leaves  the  circulating 
boost  pump  into  a  fixed  speed  centrifugal  pressure  pump  which  increases  the  pressure  of  the 
process  fluid  from  0,7  to  4.1  MPa  (100  to  595  psia).  We  always  operate  at  pressures  above 
3.5  MPa,  the  critical  pressure  of  nitrogen;  therefore,  two-phase  composition  of  the  nitrogen 
is  avoided. 
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Thermal  energy  from  three  sources  is  added  to  the  process  fluid:  (1)  the  ambient 
surroundings,  (2)  a  large  counterflow  plate-fin  heat  exchanger,  and  (3)  a  shell  and  tube  heat 
exchanger  with  steam  on  the  tube  side.  The  main  heat  exchanger  is  a  five  pass,  plate-fin 
type,  constructed  of  aluminum  and  insulated  with  a  30  cm  thickness  of  polyurethane  foam. 
The  first  pass  in  the  main  heat  exchanger  adds  heat  to  the  fluid  through  a  counterflow 
process  with  the  returning  fluid.  Downstream  of  the  first  pass  in  the  main  heat  exchanger, 
a  steam  heat  exchanger  is  used  as  a  trim  heater  to  control  the  temperature  in  the  test 
section.  The  typical  test  section  temperature  is  288  K.  The  process  fluid  returns  through 
the  second  pass  of  the  heat  exchanger,  counterflow  to  the  first  pass,  removing  heat  from  the 
return  process  stream.  Two  additional  passages  are  used  to  exchange  heat  between 
subcooled  liquid  nitrogen  and  the  returning  process  stream,  removing  additional  thermal 
energy  from  the  process  stream.  One  passageway  is  closed  and  used  as  a  gas  thermometer 
(with  pressure  sensor)  for  heat  exchanger  monitoring  and  control  purposes.  If  the  pressure 
in  the  gas  thermometer  remains  constant  during  a  test  point,  then  the  temperature  is 
constant  and  the  rate  of  heat  exchange  is  stable. 

Process  fluid  flow  is  controlled  by  the  operation  of  an  expansion  valve  downstream 
of  the  main  heat  exchanger,  in  conjunction  with  the  variable  speed  circulating  pump. 
Additional  thermal  energy  is  removed  downstream  of  the  expansion  valve  with  the  liquid 
nitrogen  bath  subcooler  mentioned  in  the  liquid  system  description.  At  this  point  in  the  flow 
loop,  the  process  fluid  is  Hquid  nitrogen  and  the  load  cell  and  reference  described  earlier 
are  used  to  complete  the  flow  loop.  To  maintain  stable  operation  of  this  system,  a  minimum 
flow  rate  of  approximately  0.5  kg/s  must  be  maintained  and  pressure  drop  in  the  gas  loop 
must  be  minimized  (<  276  kPa,  40  psia). 

This  is  a  complex  thermodynamic  cycle,  and  even  though  it  is  highly  instrumented, 
it  is  difficult  to  quantify  what  is  occurring  at  every  element  of  the  loop.  The  heat  transfer 
in  the  main  heat  exchanger  is  the  most  difficult  to  describe,  but  stability  in  the  enclosed  gas 
thermometer  reflects  a  stable  heat  exchange,  that  is,  steady  state  flow,  which  is  maintained 
throughout  the  test  point.  If  this  stabiHty  is  maintained,  the  mass  passing  through  a  meter 
in  the  gas  portion  of  the  loop  is  equivalent  to  the  mass  being  collected  in  the  weigh  tank 
during  a  data  point.  Because  we  can  not  quantify  all  uncertainties  in  this  cycle,  the  gaseous 
nitrogen  flow  loop  is  used  only  for  applied  research  and  not  for  calibrations. 

4.2   Uncertainty  in  Gaseous  Nitrogen  Mass  Flow  Measurement 

We  measure  mass  flow  in  the  gaseous  nitrogen  flow  loop  with  either  the  weigh  tank 
described  in  Section  3.1  or  by  a  gas  turbine  meter.  The  uncertainties  associated  with  each 
are  discussed. 

42.1  Weigh  Tank  Uncertainty 

The  uncertainties  in  measuring  mass  flow  rate  with  the  weigh  tank  are  explained  fully 
in  Section  3.2.  The  only  difference  between  the  analysis  in  Section  3.2  and  the  analysis 
provided  here  is  the  possible  change  in  mass  in  the  piping  volume  between  the  test  section 
(gas  versus  liquid)  and  the  weigh  tank.   If  the  system  were  totally  stable,  the  mass  in  this 
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volume  would  remain  constant  during  a  data  point. 

As  we  explained  in  the  description  of  the  gas  flow  facility  in  Section  4.1,  the  nitrogen 
in  the  piping  between  the  gaseous  nitrogen  test  section  and  the  weigh  tank  is  changed  from 
supercritical,  ambient  temperature  nitrogen  to  liquid  nitrogen.  It  is  safe  to  assume  that 
during  a  100  s  test  point,  there  may  be  some  changes  in  density  in  subvolumes  of  this  piping. 
This  would  not  be  true  for  test  points  in  which  we  measure  nitrogen  gas  flow  with  the 
turbines.  The  turbines  are  close  to  the  orifice  meter,  and  there  is  no  phase  change  in  the 
nitrogen  fluid  between  them. 

We  measure  and  record  pressures  and  temperatures  at  many  locations  between  the 
test  section  and  the  weigh  tank  to  quantify  density  changes  in  these  subvolumes.  We 
analyzed  the  changes  in  mass  between  the  gas  test  section  and  the  weigh  tank  as  a  function 
of  these  density  measurements  recorded  over  a  period  of  several  years.  We  estimate  that 
the  uncertainty  (Type  A)  in  a  nominal  mass  of  181.4  kg  due  to  the  density  variation  in  the 
piping  is  0.03%.  Table  8  shows  the  uncertainties  in  gas  measurement  when  using  the  weigh 
tank. 

Table  8.  Total  Uncertainties  (la)  for  Mass  Measurement  in  Gaseous  Nitrogen  Loop 

Using  Weigh  Tank 


Source  of  Uncertainty: 

Type  A 

TypeB 

Combined 

Load  Cell  Sensitivity 

0.050% 

0.002% 

0.050% 

Buoyancy  Correction 

0.0001% 

0.068% 

0.068% 

Mass  between  Test  Section  and  Weigh 
Tank  (gas-liquid  Nj) 

0.030% 

Total  for  Mass  Measurement 

0.058% 

0.068% 

0.089% 

Expanded  Uncertainty,  k=2 

0.179% 

4.2.2   Gas  Turbine  Meter  Uncertainty 

The  minimum  flow  rate  at  which  the  gaseous  nitrogen  flow  loop  can  operate  is 
approximately  0.5  kg/s.  We  wanted  to  test  orifice  meters  in  the  gas  loop  at  lower  flow  rates, 
so  we  added  two  gas  turbine  meters  in  series  upstream  of  the  orifice  meter.  We  bypass  part 
of  the  flow  before  it  enters  the  turbine-orifice  meter  loop  so  that  total  system  flow  is  greater 
that  0.5  kg/s.  One  of  the  turbines  is  used  as  the  primary  measurement  device,  while  the 
other  serves  as  a  constant  reference  to  the  primary  turbine. 

To  determine  the  mass  flow  rate  in  the  gas  loop,  the  calibration  curve  for  the  primary 
turbine  is  used  to  determine  the  volumetric  flow  rate  during  the  test  point.  The  density  and 
the  mass  flow  rate  are  calculated  using  the  temperature  (T17)  and  the  pressure  (P17) 
measured  at  the  turbine  meter.  The  uncertainty  in  the  mass  flow  rate  calculated  from  the 
turbine  meter  is  a  combination  of  the  uncertainties  in  the  calibration  equation,  the  counter 
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measuring  turbine  meter  pulses,  the  mass  flow  rates  used  in  the  cahbration,  and  the  density 
at  the  turbine  meter. 

4.2.2.1   Calibration  Equation 

We  calibrated  the  turbine  meters  over  the  range  of  our  gas  flow  mass-based  system: 
0.5  to  2.3  kg/s  (1.0  to  5.0  Ib/s).  The  turbine  meters  and  associated  piping  were  then  sent 
to  an  external  facihty  to  be  calibrated  over  a  flow  rate  range  of  0.11  to  1.1  kg/s  (0.25  to  2.4 
Ib/s).  This  calibration  was  done  using  air  at  pressures  and  temperatures  similar  to  those  at 
which  our  facility  operates.  The  calibration  report  included  a  stated  flow  measurement 
uncertainty  of  0.5  percent. 

The  calibration  data  sets  from  NIST  and  the  external  laboratory  were  dissimilar  and 
difficult  to  resolve.  The  data  set  taken  at  the  NIST  laboratory  contained  454  points  and 
showed  more  scatter  than  the  40  points  from  the  external  laboratory  calibration.  We 
enlisted  the  help  of  the  Statistical  Engineering  Division  at  NIST  Boulder.  They  used  a 
nonlinear  fitting  to  simultaneously  calculate  the  parameters  of  a  two-line  fit  for  volume  flow 
rate  versus  frequency. 

The  calculation  of  uncertainty  in  the  turbine  meter  calibration  equation  is  not  only 
a  function  of  the  calibration  curve  fit,  but  the  uncertainty  in  the  mass  flow  reference  system 
of  the  facility  performing  the  calibration  and  the  density  calculation.  The  turbine  calibration 
curve  fit  relies  more  heavily  on  the  data  from  the  NIST  calibration,  and  the  standard 
uncertainty  of  the  calibration  is  a  reflection  of  the  scatter  from  the  NIST  data.  For  these 
reasons,  we  chose  the  uncertainties  associated  with  measuring  mass,  time,  and  density  at  the 
NIST  lab  to  evaluate  the  uncertainty  in  the  turbine  meters. 

The  uncertainty  in  volumetric  flow  rate  measured  by  the  turbine  meter  due  to  the 
calibration  curve  fit  alone,  at  a  nominal  value  of  1.81  kg/s,  is  0.20  percent. 

42.2.2   Mass 

The  uncertainty  in  mass  measurement  was  provided  in  section  4.2.1.  We  use  only  the 
Type  B  uncertainties  here,  because  the  Type  A  uncertainties  are  reflected  in  the  calibration 
data  scatter  and  are  included  in  the  calibration  equation  uncertainty.  Type  B  uncertainty 
for  mass  measurement  is  0.068  percent. 

4.2.2.3    Frequency 

The  turbine  frequency  calculation  is  a  combination  of  the  output  of  the  universal 
counter  and  the  timer  in  the  datalogger.  A  manual  trigger  at  the  beginning  and  the  end  of 
the  measurement  interval  is  used  to  signal  measurements  in  both  instruments.  We  totalize 
the  pulses  from  the  turbine  meter,  and  no  test  point  has  fewer  than  2500  pulses.  Because 
the  uncertainty  in  the  counter  is  1  pulse,  the  maximum  uncertainty  due  to  the  counter  is  0.04 
percent.  The  Type  B  uncertainty  in  time  measurement  was  shown  in  Section  3.2.3  to  be 
0.001  percent. 
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422.4   Density 

The  calculated  fluid  density  at  the  turbine  meter  is  used  to  determine  the  mass  flow 
rate  in  the  gas  flow  loop.  Table  9  lists  the  density  uncertainties  evaluated  at  the  conditions 
noted,  along  with  the  other  contributors  to  turbine  meter  uncertainty.  As  in  Section  3.7,  the 
expanded  uncertainty  is  computed  with  and  without  the  contribution  of  MIPROPS.  The 
uncertainty  in  our  measurements  due  to  MIPROPS  would  not  be  considered  if  all  parties 
were  using  the  same  property  library. 

Table  9.   Uncertainty  (la)  in  Mass  Measurement  Determined  by  the  Turbine  Meter 


Source  of  Uncertainty 
in  Turbine  Mass  Flow 

Nominal  Values:  T= 288.7  K,  P=3.87  MPa, 
Mass  =  181.4  kg,  Time  =  100  seconds 

Type  A 

Type  B 

Calibration  Eq. 

0.20% 

Uncert.  in  mass 

0.068% 

Pulses 

0.040% 

Density:  MIPROPS 

0.10% 

(5p/(5T  p  (T17) 

0.029% 

0.057% 

dp/d?  T(P17,Baro) 

0.026% 

0.070% 

Total  (in  quad) 

0.204% 

0.156% 

Expanded  Uncertainty,  I 

:  =  2 

0.513% 

Without  MIPROPS  unce 

Ttainty,  k  =  2 

0.473% 

The  calculated  uncertainty  in  mass  measured  by  the  turbine  meter  is  more  than  twice 
that  of  mass  measured  by  the  weigh  tank  because  it  contains  the  uncertainty  of  the  weigh 
tank  as  well  as  several  other  components  (Section  4.2.2.1).  When  turbine-based  mass  flow 
rates  are  used  in  the  calculation  of  orifice  meter  discharge  coefficients,  this  additional 
uncertainty  is  not  apparent.  When  we  use  both  the  turbine  meter  and  the  weigh  tank  to 
measure  the  same  mass  flow  rate,  the  orifice  meter  discharge  coefficients  calculated  using 
the  turbine  system  calculated  mass  and  those  from  the  weigh  tank  system  measured  mass 
are  indistinguishable.  The  major  contribution  to  the  turbine  meter  uncertainty  comes  from 
the  curve  fit  to  the  calibration  equation.  The  repeatability  of  turbine-based  discharge 
coefficients  is  0.1  to  0.15  percent. 
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43  Uncertainty  in  Orifice  Meter  Discharge  Coefficient 

The  uncertainty  in  discharge  coefficient  is  calculated  using  the  propagation  of 
uncertainties  technique  described  in  Section  2.  Most  of  the  variables  in  the  discharge 
coefficient  equation  represent  values  whose  uncertainties  are  relatively  insensitive  to  flow 
rate.  This  is  not  the  case  for  the  uncertainty  associated  with  the  differential  pressure 
transducers.  The  description  of  uncertainty  in  the  differential  pressure  transducers  in 
Appendix  E  demonstrates  this  variability.  Therefore,  it  is  impossible  to  provide  a  single 
value  for  uncertainty  in  determining  discharge  coefficients. 

The  variables  in  the  equation  are  mass,  time,  density,  differential  pressure,  and  pipe 
and  orifice  diameters.  The  only  component  of  the  discharge  coefficient  uncertainty  that  has 
not  been  explained  in  detail  in  previous  sections  or  in  the  Appendices  is  the  uncertainty  in 
density.  Table  10  outlines  the  components  of  density  uncertainty  at  the  orifice  meter. 

Table  10.  Uncertainties  (la)  in  Gas  Density  at  Orifice  Meter 


Source  of  Uncertainty  in 
Orifice  Gas  Density 

Evaluated  at  3.79  MPa,  288.7  K 

Type  A 

TypeB 

MIPROPS 

0.10% 

dp/6T  p  (T16) 

0.031% 

0.057% 

dp/6?  T  (P16,Baro) 

0.039% 

0.055% 

Total  (in  quad) 

0.050% 

0.128% 

The  variables  in  the  discharge  coefficient  equation  are  defined  in  Section  2.  The  term  CY 
is  the  product  of  the  discharge  coefficient  and  the  expansion  factor  Y.  The  equation  used 
to  calculate  the  expansion  factor  is  found  in  metering  standards  and  is  associated  with  an 
uncertainty.  The  uncertainty  we  define  for  discharge  coefficient  is  actually  the  uncertainty 
in  our  calculation  of  CY.  The  discharge  coefficient  equation  and  the  method  used  to 
evaluate  the  propagation  of  uncertainty  in  the  equation  are  shown  in  the  following 
equations: 

^Y  =  (constant)m(l-^*y'^  .g. 

Fa  d^pApy^t 
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Therefore,  uncertainties  in  the  calculation  of  discharge  coefficients  are  calculated  by  solving 
the  following  equation  with  the  appropriate  uncertainties  for  all  variables.  We  combined 
the  Type  A  and  Type  B  uncertainties  prior  to  this  calculation: 
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Tables  11  and  12  provide  a  detailed  description  of  the  components  in  the  uncertainty 
calculations  of  discharge  coefficient  for  two  differential  transducers.  Table  11  lists  the 
uncertainties  when  the  weigh  tank  is  used  for  mass  measurement.  Table  12  lists 
uncertainties  when  the  turbine  is  used  for  mass  measurement. 
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5.   Summary 

Table  13  lists  what  we  consider  to  be  important  overall  uncertainties  associated  with 
this  laboratory.  These  values  are  an  estimate  of  a  95  percent  confidence  interval  using  a 
coverage  factor  of  2.  Once  again,  values  both  with  and  without  the  contribution  of  the 
density  uncertainty  due  to  MIPROPS  are  given. 

Table  13.  Uncertainty  (la)  Summary 


Flow  Measurement  Uncertainty 

Evaluated  for  181.4  kg  mass,  100  s  (k  =  2) 

Device:   Weigh  tank 

With  MIPROPS 

Without  MIPROPS 

LNj  "^^s  flow 

0.170% 

0.170% 

LNj  volume  flow 

0.530% 

0.178% 

N2  gas  mass  flow 

0.179% 

0.179% 

Device:  Turbine 

Gas  N2  mass  flow 

0.513% 

0.473% 

The  uncertainty  in  discharge  coefficient  calculation,  a  component  of  which  is  mass 
flow  rate,  is  variable  and  is  better  illustrated  by  plotting  the  relationship  between  uncertainty 
and  pressure  differential  for  the  transducers  described  in  this  report.  Figure  3  demonstrates 
this  more  complex  uncertainty  evaluation  (MIPROPS  uncertainty  included). 
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Figure  3.        Discharge  Coefficient  Uncertainty  vs.  Differential  Pressure 
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APPENDIX  A. 


LOAD      CELL      CALIBRATION     WEIGHTS      REPORTS     AND 
UNCERTAINTIES 


The  load  cell  calibration  weights  are  four  nominal  113.4  kg  (250  lb)  brass  cylinders. 
These  weights  have  been  calibrated  three  times  by  NIST  and  by  the  Colorado  Department 
of  Agriculture.  The  calibration  reports  are  included  on  the  following  pages.  A  summary 
of  the  calibration  history  is  given  in  Table  A.l,  using  the  units  of  measurement  (lb)  provided 
by  the  calibrating  agency. 

Table  A,l.   Uncertainty  (la)  in  Calibration  Weights 


Date 

Weight  #1 

Weight  #2 

Weight  #3 

Weight  #4 

Uncert. 

2-25-70 

255.34 

247.72 

252.33 

247.53 

0.01% 

10-21-71 

255.34 

247.72 

252.33 

247.53 

0.01% 

8-2-78 

255.3371 

247.7215 

252.3294 

247.5275 

0.002% 

The  1978  calibration  performed  by  the  Colorado  Department  of  Agriculture  states 
that  the  random  uncertainty  is  negligible  and  the  uncertainty  is  from  known  bias 
uncertainties,  Type  B.  The  value  of  0.002  percent  will  be  used  as  the  uncertainty  due  to  the 
weights.  The  weights  did  not  change  over  an  eight  year  period  and  we  believe  that  they 
have  remained  constant  since  that  time.  The  weights  have  been  corrected  for  local  gravity 
and  the  buoyancy  of  air. 
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COLORADO    DEPARTMENT    OF    AGRICULTURE 

406  STATE  SERVICES  BUILDING 

1525  SHERMAN  STREET 

DENVER,  COLORADO     80203 


August  2,  1978 

REPORT  OF  TEST 

for 

Four  special  purpose  test  weights 

OWNER:  National  Bureau  of  Standards,  Cryogenic  Laboratory 

Boulder,  Colorado 
SUBMITTED  BY:   J  A  Brennan 


AGRICULTURAL  COMMISSION 

Clarence  Stone,  Center 
Chairman 

William  A.  Stephens,  Gvpsum 
Vice-C  hairman 

Ben  Eastman,  Hotchkiss 
John  L.  Mailov,  Denver 
M.  C.  McCormick,  Holly 
Elton  Miller.  Fort  Lupton 
Kay  D.  Morison.  Fleming 
William  H.  Webster,  Greeley 
Kenneth  G.  Wilmore,  Denver 


Cert.  No:  ^k^^ 


The  test  weights  described  above  have  been  compared  with  the  Standards  of  the  State 
of  Colorado  and  the  following  results  have  been  determined: 


Weight  #1 
#2 
#3 


255.3371  lbs 
247.7215 
252.3294 
2.^7.5275 


Uncertainty     0.0052  lb. 


The  uncertainty  figure  is  an   expression  of  the  overall  uncertainty   using  three 
standard  deviations  as  a   limit  to  the   effect  of  random  errors  of  measurement,    the 
magnitude  of  systematic   errors  from  known   sources  being  negligible. 


F  H  Brzot).cky ,  j;!!'hief  Metrologist 

Colorad6  Depaij^ment  of  Agriculture 

Metrology   Laboratory 

3125  Wyandot   St. 

Denver,    Colorado        80211 


THESE  CERTIFlCATiCNS  ARE  TRACEABLE  TO  THE 
fJATSONAL  BUREAU  CF  STAIJDARDS. 
ALL  CERTIFICATES  ISSUED  BY  THE  COLORADO 
■^W'^T^-NT  OF  AGRiCULTURZ-METROLOGY  LAB- 
C.RhTGRY  "expire  one  YEAR  FROM  THE  DATE  BF 
ISSUANCE. 
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Date:    October  12,   1971 
X'ot:    John  Shafer,   272.55 


Subject:    Calibration  of  Brass  Cylinders 


To:  Jim  Brennan,  275*06 


U.S.  DEPARTMENT  OF  COMMERCE 
IMational  Bureau  of  Standards 

Bouldar.    Colorado      80302 


Four  brass  cylinders  were  compared  to  weights  traceable  to  the  National 
Standards  and  were  f oimd  to  be  as  follows : 


#1 

255.3^  lbs. 

#2 

2U7.72  lbs. 

#3 

252.33  lbs. 

i^ 

2i^7.53  lbs. 

The  weights  are  identified  by  numbers  stamped  at  intervals  around  the 
circumference  of  the  cylinders . 

The  measurements  were  made  March  k,   1971* 

The  data  given  above  is  correct  to  ±0.01  percent. 

Six  hours  were  charged  to  2750161 . 
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»<cD-ui      UNITED  STATES  GOVERNMEMl 

(ii-a) 

""  Memorandum 


U.S.  DEPARTMENT  OF  COMMERCE 

NATIONAL  BUREAU  OF  STANDARDS 

Boulder,  Colorado      80302 


TO  :  John  Dean 

275. 06 

FROM     :  John  F.    Shafer 

272.  24  Q*^J 

SUBJECT:  Calibration  of  Brass  Cylinders 


DATE:     February  25,    1970 
In  reply  refer  to: 


Four  brass  cylinders  were  compared  to  weights  traceable  to  the 
National  Standards  and  found  to  be  as  follows: 


#  1 

#  2 
#3 

#  4 


255.  34  lbs. 
247.72  lbs. 
252.  33  lbs. 
247.  53  lbs. 


The  weights  are  identified  by  numbers  stamped  at  intervals  arotind  the 
circunnierence  of  the  cylinders. 

The  data  given  above  is  correct  to    ±  0.  01  percent, 

7  hours  was  charged  to  Project  No.    2750360,    work  order  No.    0156. 
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APPENDIX  B. 


NITROGEN  PROPERTY  UNCERTAINTY 


Depending  on  the  location  in  the  flow  system,  we  are  measuring  subcooled  liquid 
nitrogen  or  gaseous  nitrogen  (supercritical)  flow.  The  properties  of  these  fluids  are 
calculated  using  a  computerized  package,  NIST  Pure  Fluids  Database  (MIPROPS)  [5].  The 
nitrogen  property  information  in  the  database  was  taken  from  a  report  by  B.  A.  Younglove 
[6].  The  uncertainty  expressed  by  Younglove  for  the  density  values  of  nitrogen  were  Txr, 
therefore,  the  uncertainties  that  we  list  in  our  table  are  the  stated  uncertainties  divided  by 
2.  The  density  uncertainties  are  different  for  the  two  phases  of  nitrogen  that  are  in  our 
system.  Table  B.l  lists  the  uncertainties  in  the  density.  We  also  include  the  change  in 
density  associated  with  small  changes  in  temperature  and  pressure.  These  values  are  used 
to  evaluate  density  uncertainties  in  various  regions  of  the  flow  system  using  the  Type  A  and 
Type  B  uncertainties  associated  with  temperature  and  pressure  measurements. 

Table  B.l.   Uncertainties  (Itr)  in  Nitrogen  Properties 


Property  Uncertainties 

Liquid  Nitrogen 
T=85K,P  =  621kPa 

Gaseous  Nitrogen 
T=288.7  K,  P=3.87  MPa 

MIPROPS 

Density 

0.25% 

0.10% 

(5p/dT  p 

0.703%  per  K 

0.378%  per  K 

(5p/(5P|t 

0.0005%  per  kPa 

0.025%  per  kPa 
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APPENDIX  C.  UNCERTAINTIES  IN  ELECTRONIC  INSTRUMENTATION 

Most  of  the  uncertainties  described  in  this  Appendix  are  those  stated  by  the 
manufacturers  of  the  instruments.  As  outlined  in  Technical  Note  1297  [1],  we  have  assumed 
that  the  stated  manufacturers'  uncertainties  are  based  on  a  rectangular  distribution  which 
has  a  standard  uncertainty  equal  to  the  half-width  of  the  specified  interval  divided  by  the 
square  root  of  3.  A  description  of  the  instruments  follows,  and  Table  C.l  is  a  composite  of 
the  Type  B  uncertainties  that  are  incorporated  into  other  calculations. 

Multi-function  Datalogger 

The  datalogger  is  used  for  several  purposes.  It  measures  the  output  of  the  constant- 
current  source,  it  supplies  current  to,  and  measures  the  voltage  output  from,  two 
thermometers  located  in  the  liquid  catch  tank  ullage,  and  it  measures  elapsed  time  during 
a  test  point. 

Computer  Analog  to  Digital  Converter 

The  computer  analog  to  digital  (A/D)  converter  reads  voltage  outputs  from  as  many 
as  80  analog  transmitters.  Almost  every  transmitter  in  both  the  liquid  and  gas  phase  of  the 
flow  loop  is  connected  to  this  I/O  device.  All  channels  can  be  read  at  an  approximate  rate 
of  100  times  per  second.  Because  of  the  resolution  of  this  instrument,  it  provides  the 
greatest  contribution  to  uncertainty  of  any  electronic  instrument. 

Universal  Counter 

The  universal  counter  counts  the  total  number  of  pulses  from  any  pulse-type  meter 
(turbine)  in  the  gas  flow  system. 
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Table  C.l.     Electronic  Instruments  Uncertainties  (la)  Quoted  by  Manufacturers 


Instrument 

Nominal  Value 

Type  B  Uncertainty 

Datalogger: 
A/D  Converter 

Timer 

0.1  V 

1.0  V 

10.0  V 

Variable 

0.0040% 
0.0035% 
0.0035% 

0.001  s 

Computer 
A/D  Converter^ 

Range: 

5.12-10.24  V 

2.56-5.11  V 

1.28-2.55  V 

0.0-1.27  V 

2.89  mV 

1.44  mV 

0.72  mV 

0.036  mV 

Universal  Counter 

>2500 

1  pulse 

Constant-Current  Source 

The  constant-current  source  is  used  to  supply  current  to  the  platinum  resistance 
thermometers  (PRT)  located  throughout  the  liquid  and  gas  systems.  The  calculated  current 
is  used  during  the  calibration  of  the  PRTs  and  during  temperature  measurement  with  the 
PRTs.  The  current  is  determined  by  the  datalogger  which  measures  the  voltage  across  a 
precision  resistor;  therefore,  the  uncertainty  associated  with  the  current  to  PRTs  is  actually 
a  combination  of  the  datalogger  and  resistor  uncertainties. 

We  evaluate  the  uncertainty  in  a  1  mA  signal,  the  nominal  current  supplied  to  all 
PRTs.  The  Type  B  component  of  this  uncertainty  is  a  combination  of  the  datalogger  and 
resistor  uncertainties.  The  Type  A  component  of  this  uncertainty  is  the  standard  deviation 
of  a  sample  of  readings  from  the  constant  current  source.  These  uncertainties  will  be  used 
in  the  propagation  of  uncertainties  for  PRTs  in  Appendix  F. 

Table  C.2.   Uncertainty  (la)  in  1  mA  Current  Measurement 


Current  Uncertainty 

Type  A 

Type  B 

Evaluated:  0.1  V,  100  Q 

0.000  97% 

0.0126% 

^Uncertainties    are    calculated    based    on    voltages    at    the    nominal    value    under 
consideration. 
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APPENDIX  D.  UNCERTAINTIES  OF  QUARTZ  BOURDON  GAUGES  USED  AS 

CALIBRATION  PRESSURE  STANDARD 

Quartz  Bourdon  gauges  are  highly  accurate  pressure  standards  used  to  cahbrate 
differential  and  static  pressure  transducers.  The  differential  pressure  transducers  measure 
the  pressure  drop  across  an  orifice  plate,  and  static  pressure  transducers  measure  line 
pressure.  The  interlaboratory  device  we  use  to  calibrate  the  Bourdon  gauges  is  a  dead- 
weight tester  calibrated  by  NIST.  These  calibrations  have  shown  the  Bourdon  gauges  to  be 
extremely  stable  over  time.  The  resultant  uncertainties  from  these  calibrations  are  smaller 
than  those  stated  by  the  manufacturer,  but,  for  simplicity,  we  have  decided  to  take  a  more 
conservative  route  and  use  the  uncertainties  supplied  by  the  manufacturer.  The  contribution 
of  the  uncertainty  in  these  instruments  to  overall  uncertainty  in  pressure  measurement  is 
minimal. 

The  manufacturers'  uncertainties  for  these  instruments  are  defined  as  a  percentage 
of  the  range  of  the  instrument  and  the  reading  at  a  specific  pressure.  The  output  of  the 
instrument,  read  by  the  datalogger,  is  a  voltage  which  is  proportional  to  pressure.  As 
outlined  in  Technical  Note  1297  [1],  we  have  assumed  that  the  stated  manufacturers' 
uncertainties  are  based  on  a  rectangular  distribution  which  has  a  standard  uncertainty  equal 
to  the  half-width  of  the  specified  interval  divided  by  the  square  root  of  three.  The  Type  B 
uncertainties  (Icr)  for  these  instruments  are: 

0.0017%  Range 
0.0030%  Reading 

The  contribution  of  the  uncertainties  of  these  instruments  to  the  uncertainty  in  pressure 
measurements  are  calculated  based  on  nominal  pressures  under  consideration.  These 
nominal  pressures  are  shown  in  the  tables  associated  with  each  transducer  found  in 
Appenduc  E. 
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APPENDIX  E.  UNCERTAINTIES  IN  PRESSURE  SENSING  INSTRUMENTS 

The  uncertainty  in  pressure  measurement  has  many  components.  They  include  the 
standard  error  of  the  estimate  of  the  caHbration  equation,  uncertainties  in  the  calibration 
standard  (Appendix  D),  uncertainties  in  all  electronic  equipment  (Appendbc  C),  as  well  as 
dynamic  variability  observed  during  data  points.  The  output  from  the  static  and  differential 
pressure  transducers  are  digitized  by  the  computer  A/D  converter.  The  data  acquisition 
channels  and  cables  for  the  transducers  are  the  same  during  calibration  and  actual  operating 
conditions.   The  signal  from  the  quartz  Bourdon  gauge  is  read  by  the  datalogger. 

The  dynamics  of  pipe  flow  create  variabilities  in  the  signals  from  the  pressure 
transducers  that  are  not  present  during  calibrations.  We  evaluated  the  dynamic  variability 
of  the  transducers  by  analyzing  a  sample  of  mean  output  values  for  the  both  the  static  and 
differential  transducers.  Each  value  in  the  sample  represents  the  mean  of  the  readings 
taken  during  a  measurement  (50  to  400  readings),  and  we  analyzed  at  least  57  values  for 
each  transducer. 

Because  static  or  line  pressures  vary  little  from  data  point  to  data  point,  the  standard 
errors  of  the  means  (s/Vn)  for  the  sample  of  static  pressure  values  were  averaged.  The 
average  was  used  as  the  dynamic  uncertainty.  We  divided  the  dynamic  uncertainty  by  a 
nominal  pressure  for  each  transducer  to  determine  the  fractional  uncertainty  in  pressure 
measurement  as  a  result  of  flow  dynamics. 

A  slightly  different  approach  was  used  to  evaluate  the  dynamic  variabilities  in  the 
differential  pressure  transducers.  The  standard  error  of  the  means  of  the  differential 
pressure  measurements  were  fitted  versus  the  mean  differential  pressures.  Using  this  fit,  we 
determined  the  dynamic  uncertainty  at  specific  differential  pressures  shown  in  Tables  E.6 
and  E.7. 

Atmospheric  and  Static  Pressure 

The  barometer  which  measures  ambient  air  pressure  is  also  calibrated  with  the  quartz 
Bourdon-tube  gauge.  The  barometer  reading  is  used  for  calibrating  the  static  pressure 
transducers  and  its  uncertainty  contributes  to  the  uncertainty  in  the  static  pressure 
measurement.  The  pressure  range  measured  by  static  pressure  transducers  varies  from  3.72 
to  4.0  MPa  (540  to  580  psia).  The  uncertainty  for  these  transducers  will  be  evaluated  at  a 
nominal  pressure  of  3.87  MPa  (562  psia).  The  following  tables  outline  the  uncertainties 
associated  with  the  static  pressure  transducers. 
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Table  E.l.  Barometer  Uncertainties  (la) 


Barometer  Uncertainties 
(Evaluated  @  87.9  kPa) 

Type  A 

Type  B 

Calibration  Equation 

0.0062% 

Quartz  Bourdon  Gauge 

0.0047% 

Datalogger  A/D 

0.0014% 

Computer  A/D 

0.046% 

Total  (in  quadrature) 

0.0062% 

0.046% 

Table  E.2.   Uncertainties  (la)  in  Orifice  Meter  Static  Pressure  Transducer 


Orifice  Static  Pressure 
Uncertainties 
(P=3.87  MPa) 

Type  A 

TypeB 

Calibration  Equation 

0.0405% 

Quartz  Bourdon  Gauge 

0.0046% 

Datalogger  A/D 

0.0035% 

Computer  A/D 

0.034% 

Barometer  (87.9  kPa) 

0.0062% 

0.046% 

Flow  Dynamics 

0.002% 

Total  (in  quadrature) 

0.041% 

0.057% 
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Table  E3.  Uncertainties  (la)  in  Tarbine  Meter  Static  Pressure  Transducer 


Turbine  Static  Pressure 
Uncertainties 
(P=3.87  MPa) 

Type  A 

TypeB                    1 

Calibration  Equation 

0.0265% 

Quartz  Bourdon  Gauge 

0.0046%                     1 

Datalogger  A/D 

0.0035% 

Computer  A/D 

0.056% 

Barometer  (87.9  kPa) 

0.0062% 

0.046% 

Flow  Dynamics 

0.002% 

Total  (in  quadrature) 

0.027% 

0.073% 

Table  E.4.  Uncertainties  (la)  in  Liquid  Test  Section  Static  Pressure  Transducer 

P7  LN2  Static  Pressure 
Uncertainties  (P=586  kPa) 

Type  A 

TypeB 

Calibration  Equation 

0.061% 

Quartz  Bourdon  Gauge 

0.014% 

Datalogger  A/D 

0.0035% 

Computer  A/D 

0.034% 

Barometer  (87.9  kPa) 

0.0062% 

0.046% 

Flow  Dynamics 

0.03% 

Total  (in  quadrature) 

0.068% 

0.059% 

We  did  not  perform  the  same  dynamic  analysis  for  transducer  P9,  located  in  the  weigh  tank, 
because  we  do  not  keep  the  same  level  of  archival  data  for  P9.  However,  we  estimate  the 
results  would  be  comparable  to  those  in  P7  because  the  transducers  are  similar  and  are 
located  in  similar  environments. 
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Table  E.5.   Uncertainties  (la)  in  Weigh  Tank  Ullage  Static  Pressure  Transducer 


P9  Ullage  Static  Pressure 
Uncertainties  (P=517  kPa) 

Type  A 

Type  B 

Calibration  Equation 

0.047% 

Quartz  Bourdon  Gauge 

0.015% 

Datalogger  A/D 

0.0035% 

Computer  A/D 

0.038% 

Barometer  (87.9  kPa) 

0.0062% 

0.046% 

Flow  Dynamics 

0.030% 

Total  (in  quadrature) 

0.056% 

0.062% 

Differential  Pressure  Transducers 

Differential  pressure  transducers  measure  the  pressure  drop  across  an  orifice  plate 
and  are  designed  for  certain  pressure  ranges.  The  pressure  range  required  is  a  function  of 
pipe  size,  beta  ratio,  and  flow  rate.  From  one  to  four  differential  pressure  transducers  are 
used  at  one  time  and  are  read  simultaneously.  We  evaluated  the  uncertainty  of  several 
transducers  that  cover  various  differential  pressure  ranges  and  have  selected  two  transducers 
with  the  greatest  uncertainty  to  provide  a  conservative  estimate  of  the  uncertainty  of  all 
transducers  of  similar  use  and  range.  The  uncertainty  of  the  transducers  also  depends  upon 
the  differential  pressure  at  which  it  is  operating.  We  evaluated  each  transducer  at  three 
differential  pressures. 

Tables  E.6  and  E.7  list  the  uncertainties  for  the  differential  pressure  transducers. 
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APPENDIX  F.  UNCERTAINTIES  IN  TEMPERATURE  MEASUREMENT 

We  use  platinum  resistance  thermometers  (PRTs)  to  measure  temperatures  critical 
to  accurate  flow  measurement.  They  are  calibrated  against  our  interlaboratory  standard 
platinum  resistance  thermometer  (SN  480)  that  was  calibrated  by  the  NIST  Thermometry 
Group,  Process  Measurement  Division.  The  calibration  of  system  thermometers  consists  of 
supplying  a  constant  current  to  the  PRTs  in  series  with  SN  480  and  the  precision  resistor 
(for  the  current  measurement).  The  PRTs  are  immersed  in  baths  at  various  temperatures 
which  reflect  the  temperatures  at  which  the  PRTs  will  operate  within  the  flow  system.  The 
voltage  across  the  PRTs  is  amplified  before  being  digitized  by  the  computer's  A/D 
converter.  The  voltage  across  SN  480  is  digitized  by  a  high  resolution  multi-function 
datalogger  which  transmits  the  digitized  voltage  to  the  computer  via  an  IEEE  488  databus. 
The  voltage  across  the  precision  resistor  used  for  calculating  the  current  is  digitized  by  the 
datalogger,  also. 

A  Unear  least-squares  regression  fit  of  the  voltage  from  the  PRTs  and  the 
temperature  calculated  from  SN  480  is  performed  over  the  temperature  range  of  interest. 
This  calibration  is  used  in  the  data  reduction  of  temperature  measurements.  In  the  flow 
measurement  process,  the  thermometers  remain  in  the  same  series  and  on  the  same  A/D 
converter  channels  as  used  for  calibration.  SN  480  is  the  only  item  removed  from  the 
system.  The  uncertainty  associated  with  the  thermometers  includes  the  uncertainty  due  to 
the  calibration  equation  as  well  as  the  propagation  of  uncertainties  of  current  and  voltage 
in  the  calibration  equation,  and  the  uncertainty  of  SN  480.  The  equations  used  to  perform 
the  propagation  of  uncertainty  for  the  thermometers  are  listed  below: 

T  =  A  +  BiV/I)  ,  (Fl) 


11  =  :?,  (F2) 

bV       I 


,2    _ 


uW  =  B 


11  =  ^^,  (F3) 

6/  /2 


(F4) 


We  identify  six  PRTs  that  we  consider  important  to  our  determination  of  system 
uncertainty.  Any  uncertainties  in  temperature  measurement  ultimately  contribute  to 
uncertainties  in  density  calculations.  For  this  reason,  temperature  uncertainties  have  been 
evaluated  in  kelvins.  These  values  are  then  used  to  evaluate  a  component  of  density 
uncertainty. 
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We  consider  three  of  these  thermometers  critical  to  accurate  flow  measurement: 
T16,  the  orifice  meter  thermometer,  T17,  the  turbine  meter  thermometer,  and  TIO,  the 
liquid  test  section  thermometer.  TIO  is  calibrated  between  76  K  and  85  K,  and  T16  and  T17 
are  calibrated  between  273.15  K  and  305  K.  Tables  F.l  through  F.3  illustrate  all  known 
components  of  uncertainty  for  these  thermometers. 

Table  F.l.  Uncertainties  (la)  in  Orifice  Meter  PRT,  T16 


T16-Orifice  PRT 
(Evaluated  at  288.7  K) 

Type  A 

TypeB 

Standard  PRT  (SN  480) 

0.0001% 

Voltmeter  for  SN  480 

0.0014% 

Calibration 

0.028% 

Propagation  of  Uncertainties 
(T=A+B(V/I)) 

0.0037% 

0.0528% 

At  288.7  K 

0.082  K 

0.152  K 

Table  F.2.  Uncertainties  (la)  in  Turbine  Meter  PRT,  T17 

T17-Turbine  PRT 
(Evaluated  at  288.7  K) 

Type  A 

Type  B 

Standard  PRT  (SN  480) 

0.0001% 

Voltmeter  for  SN  480 

0.0014% 

Calibration 

0.028% 

Propagation  of  Uncertainties 
(T=A+B(V/I)) 

0.0037% 

0.0527% 

At  288.7  K 

0.082  K 

0.152  K 
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Table  F.3.   Uncertainties  (la)  in  Liquid  Test  Section  PRT,  TIO 


TlO-Liquid  Test  Section  PRT 
(Evaluated  at  85  K) 

Type  A 

Type  B 

Standard  PRT  (SN  480) 

0.0001% 

Voltmeter  for  SN  480 

0.0014% 

Calibration 

0.041% 

Propagation  of  Uncertainties 
(T=A+B(V/I)) 

0.008% 

0.021% 

At85K 

0.036  K 

0.018  K 

Three  other  thermometers  provide  secondary  contributions  to  the  buoyancy 
corrections  used  in  the  weigh  tank  system.  T13  measures  the  temperature  of  the  liquid 
collected  in  the  weigh  tank.  The  weigh  tank  temperature  is  used  in  the  calculation  of 
density  and,  thereby,  volume  of  the  collected  liquid.  TS71  and  TS72  measure  the 
temperature  in  the  ullage  gas  surrounding  the  weigh  tank  and  are  used  in  the  calculation 
of  ullage  gas  density.  The  buoyancy  of  the  liquid  in  the  ullage  gas  is  the  product  of  the 
Hquid  volume  and  the  ullage  gas  density.  Though  accurate  measurement  of  the  ullage  gas 
temperature  is  important,  its  contribution  to  the  total  system  uncertainty  is  secondary. 

Thermometers  T13,  TS71,  and  TS72  are  located  in  the  catch  tank  and  are  not  as 
readily  accessible  as  those  located  in  the  flow  loop.  For  that  reason,  they  are  not  calibrated 
as  frequently.  The  residual  standard  deviation  of  the  original  calibration  of  T13  was  0.4  K. 
Because  of  the  long-term  stability  of  PRTs  and  the  fact  that  T13  remains  in  a  protected 
environment  free  from  mechanical  shock,  we  think  that  the  performance  has  not  changed 
significantly  since  the  original  calibration.  We  chose,  however,  to  use  a  slightly  greater  value 
of  0.5  K  as  a  conservative  Type  B  estimate  of  the  total  uncertainty  in  this  temperature 
measurement. 

The  manufacturer  of  thermometers  TS71  and  TS72  provided  a  one-point  calibration 
by  measuring  the  resistance  at  273.15  K.  Using  the  manufacturer's  specifications  for 
platinum  resistance  thermometers  of  this  type,  we  derived  a  general  second  order  calibration 
curve  whose  independent  variable  is  the  ratio  of  the  resistance  measured  at  operating 
temperatures  to  the  resistance  at  273.15  K.  A  test  of  this  regression  fit  at  liquid  nitrogen 
temperatures  indicated  a  possible  0.5  to  0.8  K  error  with  regards  to  temperatures  measured 
by  SN  480.  We  will  use  a  value  of  1  K  as  a  conservative  Type  B  estimate  of  the  total 
uncertainty  in  these  temperature  measurements.  Table  F.4  lists  the  uncertainties  associated 
with  these  thermometers. 
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Table  F.4.  Uncertainties  (la)  in  Weigh  Tank  and  Ullage  PRTs 


Secondary  PRT 
Uncertainties 

Type  A 

TypeB 

T13 

0.5  K 

TS71,  TS72 

1.0  K 
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APPENDIX  G. 


UNCERTAINTIES  IN  DIMENSIONAL  MEASUREMENTS 


The  uncertainties  in  measurement  of  pipe  diameter  and  orifice  diameter  are  listed 
in  Table  G.l.  These  values  will  be  used  in  the  propagation  of  uncertainty  for  discharge 
coefficients.  The  magnitudes  of  the  measurement  uncertainties  are  shown  in  parenthesis 
and  are  evaluated  at  nominal  values  of  10.2  cm  (4  in)  for  pipe  diameter  and  6.35  cm  (2.5 
in)  for  orifice  diameter. 

Table  G.l.   Dimensional  Uncertainties 


Dimensional  Uncertainty 

TypeB 

Pipe  Diameter  (0.0254  mm) 

0.025% 

Orifice  Diameter  (0.002  54  mm) 

0.004% 

*U.S.  GOVERNMENT  PRINTING  OFFICE:1994-573-013/00021 
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